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£T) , Abstract 

Q^ ' The low-energy parameters describing the np scattering in the s S± + 3 Di 

partial waves, the deuteron parameters, and their relations are discussed. 
, These parameters can be determined quite accurately in the energy-dependent 

Nijmegen partial- wave analyses of the np data and are also given by the high- 
quality Nijmegen potentials. The newest values for these parameters are pre- 
sented. 

1 Introduction 

The energy-dependent partial-wave analyses (PWA) of the NN scattering data as 
performed in Nijmegen produce very good descriptions of these NN data below 
T L = 350 MeV. We will consider explicitly the Nijm PWA93, which has x 2 /N d = 0.99 
with respect to the scattering data The deuteron binding energy B = 2.224575(9) 
MeV was included in these analyses as an experimental datum [0]. Because we have 
an energy-dependent solution of the scattering, we can produce accurate values for 
the low-energy scattering parameters and for the deuteron parameters. The PWA 
produces statistical errors on these parameters. 

In Nijmegen we have also constructed high-quality (HQ) NN potentials 0. These 
HQ potentials give a very good fit to the NN scattering data, including the deuteron 
binding energy. For each of these potentials we find in a direct comparison with the 
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NN scattering data that x 2 /N& = 1.03. Because the fit with these potentials is so 
good, we may speak here of alternative PWA's. Comparing the values of the low- 
energy parameters and the deuteron parameters as obtained in the various analyses 
gives us an estimate of the systematic error. 

We can use these HQ potentials also to compute deuteron parameters that cannot 
be determined from the NN scattering data, such as the d-state probability p d , the 
electric quadrupole moment Q, and the mean-square deuteron radius (r 2 ) 1//2 . It turns 
out that also these parameters seem to be almost uniquely determined. 

2 Scattering Matrix 

The wave function is a solution of the relativistic Schrodinger equation 

(A + k 2 )ip = 2mVip , 

with m the reduced mass. The relativistic relation between the cm-momentum k and 
the cm-energy E is 

E = ^k* + Ml + y/k 2 + Ml - (M p + M n ) , 

where M p and M n are the proton and neutron mass. The asymptotic part of the 
wave function defines the partial-wave S matrix, which is related to the K matrix by 

S= (l + iK)/(l-iK) . 

For the spin-singlet waves ( 1 So, 1 P 1 , . . .) and the uncoupled spin-triplet waves 
( 3 Po, 3 -Pi> 3 D 2 ,...), there is a simple relation between the phase shift 5 and the 
5" and K matrices, viz. 

K = taji5 and S = e 2iS . 

One can show that for low momenta 5g ~ k 2£+l . 

For the triplet coupled waves ( 3 Si + 3 Di, 3 P2 + 3 F 2 , . . .), there are two different 
parametrizations in use, the Stapp parametrization [|J] and the Blatt and Biedenharn 
parametrization [H. 

Stapp parametrization [|J 

The (2 x 2) S matrix is written in terms of the nuclear-bar phase shifts 5 and e as 

/ e^- 1 • \ / cos2£j ism2ej \ ( e^- 1 • \ 
\ ■ e iSj+1 ) \ i sin 2e j cos 2ej J \ ■ e iSj+1 ) 
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The advantage of this parametrization is that the behavior of the nuclear-bar mixing 
angle e for low momenta is fine, in contradistinction with the behavior of the Blatt 
and Biedenharn || mixing angle e (see below). For low momenta 

r 1,2,7-1 r u 2J + 3 ? t.2J+l 

Below the threshold E = the analyticity properties of the nuclear-bar phases are 
not so nice. 



Eigenphases (Blatt and Biedenharn) [BJ 

The (2 x 2) 5 matrix, and therefore also the K matrix, can be diagonalized by a real 
orthogonal matrix 

q_( cose j — sin e j 
\ sine j cose j 

with 

S = S diag O- 1 and K = O K diag O' 1 . 

The eigenvalues define then the eigenphases <5j_i and Sj + i. The difficulty with this 
parametrization is that for low momenta 

5j_ x ~ k 2J ' 1 , 5 J+1 ~ A; 2J+3 , but ej~k 2 . 

For J > 1 and small values of k the mixing parameter e,j can have already a size- 
able value, when the phases are both still practically zero. On the other hand, the 
analyticity properties of the eigenphases are nicer than of the nuclear-bar phases. 




3 Analyticity 

The S'-matrix elements are analytic functions of the complex energy E or of the 
complex momenta k 2 . They have the so-called unitarity cut from E = to oo 
along the positive real E axis. There are additional right-hand cuts starting at 
the pion-pro duct ion threshold E cm ~ 135 MeV and left-hand cuts due to meson 
exchange. The OPE cut starts at E = -ml /(AM) = -5 MeV, the TPE cut starts at 
E = —m 2 /M = —20 MeV, where M = the nucleon mass and = the pion mass. 
When there is a bound state present the S'-matrix elements have a pole at the position 
of the bound state. In the 3 Si+ 3 Di coupled channel there is therefore the deuteron 
pole at E = -2.224575(9) MeV. 

It is useful to define the effective-range function || [7j 

M(k 2 ) = k'^ITW . 

This effective-range function is again an analytic function of E or k 2 with right- and 
left-hand cuts as in the S'-matrix elements. However, the unitarity cut is removed 
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and there is no pole at the position of the deuteron. The effective-range function 
M(k 2 ) is therefore regular in a circle of radius 5 MeV around the origin, because then 
the OPE left-hand cut is reached. This effective-range function can be expanded in a 
power series around the origin with a radius of convergence of 5 MeV. One gets then 
the effective-range expansion 

M = M + M x k 2 + M 2 k 4 + M 3 k 6 + ... 

For practical purposes the radius of convergence of this effective-range expansion is 
pretty small (only 5 MeV). It is possible to define modified effective-range functions 
in which some of the left-hand cuts are removed, such as the OPE cut. The radius 
of convergence is then 20 MeV. 

In the Nijmegen PWA we do not use the modified effective-range expansions || |9| 
anymore, but our methods, however, are quite similar. 



4 The Deuteron in Scattering 

The deuteron appears as a pole in the partial-wave S matrix for the coupled 3 Si+ 3 Di 
partial waves at k — ia. The binding energy is given by 



B = M p + M n — \JM 2 -a 2 - sjMl - a 2 = 2.224575(9) MeV . 

The radius R of the deuteron is then 

R=l/a = 4.318946 fm . 

Here we use relativistic kinematics. The value of R in the case of non-relativistic 
kinematics is significantly different (R n . v . = 4.317667 fm). 

The partial- wave S matrix we write in the Blatt and Biedenharn parametrization 

S -°( S ° si) ' 1 - 

The matrix elements of the matrix O are well behaved in the neighborhood of the 
deuteron pole. This allows us to define at the pole 

7] = - tan ei = 0.02543(7) . 

The values quoted are from the Nijmegen 1993 PWA [0. One of the eigenvalues (we 
take So) has a pole at k = ia. We may write therefore 

N 2 

S = i — r + (regular function of k) . 

a + ik 

This expression is valid in the neighborhood of the deuteron pole. The residue at 
the pole is 

N 2 = 0.7830(7) far 1 . 
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It is very convenient to define the effective range pd = p(—B, —B) at the deuteron 
pole by 

p d = R- 2/N 2 = 1.765(2) fm . 

Only 3 independent deuteron parameters B, rj, and pd can be determined from the 
np scattering data. 



5 The Deuteron as a Bound State 

The deuteron is a bound state in the coupled 3 Si + 3 Di two-nucleon system. Let us 
neglect the possible existence of other coupled channels, such as A A, Q 6 , etc. In that 
case the deuteron wave function is 

, _<r) m w(r) m 
r — J'oii H — J^ll ) 

where 3^LSj m are the simultaneous eigenfunctions of the operators L 2 , S 2 , J 2 , and 
J z . We assume the wave function to be properly normalized 

POO 

/ dr (u 2 + w 2 ) = 1 . 
J o 

The asymptotic behavior of the wave functions for large values of r are 

u{r)~A s e- ar and w(r) ~ A D e~ ar {l + 3 (R/r) + 3 (R/r) 2 } . 

The asymptotic normalizations are 

A s = 0.8841 fm~ 1/2 and A D = 0.0224 fm~ 1/2 , 

with Nj = A 2 S + A 2 D = 0.7821 fm -1 . The d/s ratio r, = A d /A s = 0.02534. The values 
quoted are from the Nijm I deuteron ||. 

For the ordinary energy-independent local potentials and for velocity-dependent 
potentials one can show that Nd = N p , i.e., the normalization as defined by the 
wave functions is equal to the normalization as defined by the residue |10 |. In that 
case there is a trivial relation between the pole parameters determined in scattering 
and the asymptotic wave function. When we have energy-dependent potentials with 
dll/dk 2 0, then one can show that 

poo / qtj \ 
(N d /N p ) 2 = 1 - j o drUd^M ■ 



5 



6 Effective- Range Functions 



The (2 x 2) matrix effective-range function is M(k 2 ) = k +1 K _1 k . This matrix 
effective-range function can be expanded in a power series around the origin, but of 
course also around any other point in its analyticity domain, for example around the 
deuteron energy. 

One can do the expansion around k 2 = 0, then 

M(k 2 ) = M + Mik 2 + M 2 k A + M 3 k 6 + M 4 k 8 + . . . , 

but one could do this also around k 2 = —a 2 , then 

M(k 2 ) =m + m x (k 2 + a 2 ) + m 2 (k 2 + a 2 ) 2 + ... . 

This means that 

m = M - M x a 2 + M 2 a 4 - M 3 a 6 + M 4 a 8 

It is surprising, but also alarming, to discover that in order to get mo with sufficient 
accuracy one needs to include the term Mqo. 12 in the right-hand side. 

When we use the Blatt and Biedenharn parametrization in eigenphase shifts and 
coupling parameters, then we can define the effective-range functions 

Fj.^k 2 ) = k 2J ' x cot 5j-i , F J+l (k 2 ) = k 2J+3 cot5 J+1 , and F £ (k 2 ) = 2k 2 cot2e,j . 

These functions have the nice analyticity properties that make series expansions 
around the origin (k 2 = 0) or around the deuteron (k 2 = —a 2 ) possible. 

Using the Stapp parametrization we could define the effective-range functions 

= 2k 2J - x cot 25j_i , Fj+i = 2k 2J+3 cot 25 J+1 , and F £ = 2k 2J+1 cot 2ej . 

Again, the analyticity properties of these functions are not so nice. 

7 Effective Ranges 

It is customary to define various effective ranges such as the standard effective range 
r = p(0, 0), the mixed effective range p m = p(0, —B), and the deuteron effective range 
Pd = p(—B, —B). Let us look at the effective-range function F (k 2 ) = kcotS for the 
eigenphase 5q, also called the 3 5*i phase shift. Expanding in a power series around 
the origin gives 

F (k 2 ) = -1/a + ^rk 2 + v 2 k A + v 3 k e + . . . 

where a = scattering length. 

Some of the low-energy scattering parameters are given in Table |l| Remarkable is 
the agreement even in the parameters v 2 to v§. Numerical difficulties in our calcula- 
tions prevented us to calculate accurately vq. We think that the values for f 5 are still 
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accurate, but we are not sure. That there is such an enormous agreement in these 
expansion coefficients comes in our opinion from our correct treatment of the OPE. 
It is this OPE that determines the fast energy dependence of the scattering matrix. 
The heavier-meson exchanges give rise to a much slower energy dependence. 

It is easy to see that, because S = (1 + iK )/(l — iK ), one obtains a pole 
at k 2 = —a 2 in the S-matrix elements when K (— a 2 ) = —%. This implies that 
F (-a 2 ) — kK^ 1 — (iat)(i) — —l/R. 
A series expansion around the deuteron is 

F (k 2 ) = -l/R + \p d (k 2 + a 2 ) + w 2 (k 2 + a 2 ) 2 + ... 

The mixed effective range p m = p(0, —B) is defined by 

-1/R = -l/a + \p m (-a 2 ) . 

Then 

p m = 2R(1 - R/a) = 1.754(2) (3) fm . 

The first entry is the statistical error, the second one is the systematic error. 
We can make the expansions: 

r = Pm + 2a 2 f 2 — 2a A v 3 + 2a 6 v 4 — . . . 
Pd = p m - 2a 2 v 2 + Aa 4 v 3 - 60;% - . . . 

The terms with v 6 in the above expansions are still significant. It turns out that 

r p m — 0.001 fm and Pd ^ Pm + 0.010 fm . 

Using the value for p m gives 

r = 1.753 fm and p d = 1.764 fm . 

At this point it becomes important to see how many independent low-energy param- 
eters and deuteron parameters there are. 



Table 1: Low-energy scattering parameters for the 3 5i eigenphase. 





a 


r 


v 2 


^3 


v 4 


^5 


PWA 


5.420(1) 


1.753(2) 


0.040 


0.672 


-3.96 


27.1 


Nijm I 


5.418 


1.751 


0.046 


0.675 


-3.97 


27.2 


Nijm II 


5.420 


1.753 


0.045 


0.673 


-3.95 


27.0 


Reid93 


5.422 


1.755 


0.033 


0.671 


-3.90 


26.7 
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Table 2: Deuteron parameters. 





V 


Pd 


Pd 


(r 2 ) 1 / 2 


Qo 


PWA93 


0.02543 


1.7647 








Nijm I 


0.02534 


1.7619 


5.664 


1.9666 


0.2719 


Nijm II 


0.02521 


1.7642 


5.635 


1.9675 


0.2707 


Reid93 


0.02514 


1.7688 


5.699 


1.9686 


0.2703 



The binding energy B is directly related to the deuteron radius R. From the 
low-energy scattering follows the scattering length a. Using R and a gives the mixed 
effective range p m . This determines then the values of the effective range r and the 
deuteron effective range pd- This in turn determines the residue iV 2 at the pole. 

The d/s ratio r\ is also an independent quantity. We find therefore only 3 inde- 
pendent quantities: B, r/, and a. 

8 Multi-energy Solutions for S 

As already pointed out in the introduction we have 4 excellent multi-energy 
parametrizations for the S matrix. The first one is our energy-dependent phase-shift 
analysis Nijm PWA93. We feel that the parameters determined in this analysis [l[ are 
probably the most accurate. Then there exist 3 HQ potentials [0] (Nijm I, Nijm II, 
and Reid93) which can be considered as alternative PWA's with a slightly higher 
Xmm than Nijm PWA93. In Table |TJ we presented already some of the effective-range 
parameters. In Table ^| we present the various deuteron parameters. The deuteron 
parameters rj and pd can be determined in PWA's as well as from the HQ poten- 
tials. The d-state probability pd, the mean-square radius (r 2 ) 1 / 2 , and the electric 
quadrupole moment can only be determined for the HQ potentials. 

9 Conclusions 

In this contribution we presented the various low-energy parameters and the deuteron 
parameters. We would like to stress once more that effective-range expansions are not 
very accurate, unless one takes many terms into account. A similar situation exists 
in the 1 5'o channel. Also there the effective-range expansion is not accurate enough 
to describe the data with sufficient accuracy ||. As a result, the effective-range 
parameters have only a very limited usefulness. The following low-energy parameters 
and deuteron parameters we would like to recommend. The errors are educated 
guesses for the total error, including statistical as well as systematic errors. 
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scattering length a = 5.4194(20) fm. 
effective range r = 1.7536(25) fm. 

The values of a and r are strongly correlated due to the accurate value of the 
deuteron binding energy. A large value of a implies a large value of r. 

shape parameters t> 2 = 0.040(7) fm 3 , v 3 = 0.673(2) fm 5 , t> 4 = 3.95(5) fm 7 , and 
v 5 = 27.0(3) fm 9 . 

deuteron binding energy B = 2.224575(9) MeV. 
deuteron radius (rel) R = 4.318946 fm. 
d/s-ratio rj = 0.0253(2). 

This value is in good agreement with the recent determination by Rodning and 



Knutson |TT| of r] = 0.0256(4). 

deuteron effective range pd = 1.765(4) fm. 

residue A 2 = 0.7830(15) fm. 

asymptotic normalizations A s = 0.8845(8) fm~ 1/2 and A D = 0.0223(2) fur 1/2 . 

d-state probability pd = 5.67(4)%. 

mean- square radius (r 2 ) 1/2 = 1-9676(10) fm. 

This must be compared with the recent experimental determination in atomic 
physics [IT 7 !], which results in: (r 2 ) 1//2 = 1.971(6) fm. For good discussions see 
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• electric quadrupole moment Q = 0.271(1) fm 2 . 

This value disagrees with the experimental value Q = 0.2859(3) fm 2 . The 
difference Q — Qq ~ 0.015 fm 2 needs to be explained in terms of meson exchange 
currents, relativistic effects, contributions of Q 6 states, etc. 

We would like to point out here the existence of the NN-OnLine facility of the 
Nijmegen group [|l(|, which is accessible via the World-Wide Web. There one can 
obtain the various Nijmegen e-prints, the fortran codes for some of the Nijmegen 
potentials, the deuteron parameters and the deuteron wave functions, the phases ob- 
tained from the Nijmegen PWA and from the Nijmegen potentials, and predictions 
for many of the experimental quantities. A direct comparison of these predictions 
with the Nijmegen NN-data base is also possible. 
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